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Starting from the Kubo formula we calculate the shear 
viscosity in hot 0* theory nonperturbatively by resumming 
ladders with a real-time version of the Bethe-Salpeter equa- 
tion at finite temperature. In the weak coupling limit, the 
generalized Fluctuation-Dissipation Theorem is shown to de- 
couple the Bethe-Salpeter equations for the different real-time 
components of the 4-point function. The resulting scalar in- 
tegral equation is identical with the one obtained by Jeon 
using diagrammatic "cutting rules" in the Imaginary Time 
Formalism. 

A systematic field theoretical calculation of the vis- 
cosity from the Kubo formula involving the stress tensor 
correlation function |l|] has been the subject of intense 
theoretical interest ffla]. First results from a one- loop 
calculation [0-B of the viscosity in hot 0^ theory were 
incomplete because an infinite number of higher order 
diagrams in the loop expansion contribute to the lowest 
order in powers of the coupling constant [p|. Their re- 
summation requires a nonperturbative calculation. Jeon 
[pj identified the dominant class of diagrams contribut- 
ing to the viscosity at leading order in the coupling con- 
stant and resummed them by using so-called diagram- 
matic "cutting" rules in the Imaginary Time Formalism 
(ITF). In this paper we reproduce his results using an al- 
ternative and, we believe, more economic approach based 
on a solution of the Bethe-Salpeter equation in real time 
using the Closed Time Path (CTP) formalism B-§. Our 

I calculation draws heavily on recently derived [3| general 

I relations among the different thermal components of real- 

|time n-point Green functions. 

The CTP formalism has the advantage over the ITF 

' that it is easily generalized to non-equilibrium situations 
and that it does not require analytic continuation of the 
n-point Green functions to real time at the end of the 
calculation. It's popularity is, however, decreased by the 
doubling of degrees of freedom and the resulting compli- 
cated matrix structure of the Green functions. We begin 
by reviewing some important relations among the compo- 
nents of the n-point Green functions in the CTP forma- 
lism which help to reduce the complexity of the real-time 
calculation before evaluating any Feynman diagrams. In 
the so-called single-time representation of the CTP for- 
malism the n-point Green function is defined as 

Ga,...a„ (1, . . . , n) = (-z)«-i(rp[0,, (1) • • • 0,„ (n)]) . (1) 

Here (. . .) denotes the thermal expectation value, the 
numbers l,...,n stand for Minkowski space coordi- 
nates xi, . . . ,Xn, Tp represents the time ordering op- 



erator along the closed time path and corresponds to 
normal (antichronological) time ordering of operators 
with time arguments on its upper (lower) branch, and 
ai,a2,...,a„ — 1,2 indicate on which of the two 
branches the fields are located. (For details see, e.g., 
[||j|].) Using the KMS condition ^^ one derives in mo- 
mentum space B 



r, 
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(-l)"-lGa,a....a„(fcl,fe,...,fc„) J] 
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where ki+k2+ ■ ■ ■ +kn = 0, the star denotes complex 
conjugation, /3 is the inverse temperature, and Hi — 2,1 
for Ui = 1,2, respectively. The "physical" or r/a repre- 
sentation is defined by setting ||] 

(l)a{x) = (t>i{x)-(l)2{x), (t)r{x) = , (3) 



and writing the n-point Green function as 

^ai ...ctji l^-L, . . . , n 



nn^ — 1 

-—^{Tp[4>c.^{l)---4>c.Jn)]). (4) 



Here ai, . . . ,an — a,r, and Ur is the number of r indices 
among (ai, 012, ... , a„). One can show M that 



GQa...a(l, •■•,") = 0, 



(5) 



that the functions with only one r index, Gra...a, Gara...a, 
• • • I Ga...ar, E^re the fully retarded n-point Green func- 
tions defined in O], and that (except for Grr...r) all other 
components G„i...a„(l, . . . , n) involve both retarded and 
advanced relations among their n time arguments. The 
simplest example is the 2-point function: 



Aaa(fc) - 0, A™(fc) = A'-^*(fc), Aar{k) - 

A,,(fc) = (l+2n(fc"))(A™(fc) - Aar{k)) 



1 adv 



(fc), (6) 
(7) 



Eq. (M) is the fluctuation-dissipation theorem |jT3], with 
n(fc°) = l/{e^'^ ~1)- The retarded and advanced 2-point 
functions satisfy 

A:„(fc) = A„(fc), (8) 

A™(-fc) = Aar{k), Arr{-k) = A^^(fc). (9) 

The the single-time and r/a representations of the n- 
point Green functions are related by H 
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where repeated indices are summed over and 

Qal = —Qa2 — Qrl — Qr2 — —7= 



(11) 



are the four elements of the orthogonal Keldysh trans- 
formation for 2-point functions |1C[] . n-point Green func- 
tions involving only a single field are symmetric under 
particle exchange: 



'-^ ...ai...aj ...\- ■ ■ 5 "'«; ■ ■ ■ i '^ji • • • ) 
^ ...otj ...ai...\- ■ ■ 5 "-J 5 ■ ■ ■ 7 '^i: • • • ) 



(12) 



In the following we study massless (f>'^ theory in the 
weak coupling limit: 



A = ^(9M<^)'+^<^^ A«i. 

The Kubo formula for the shear viscosity [0,|| 



(13) 



77 =£ lim d^xe'P-^7:i„,{x)n'"'m (14) 

involves the correlation function of the traceless pressure 
tensor Trjm which in the thermal rest frame reads 



TTlmix) = {SliSmj - ^SlmSij) di(t)dj4> 



(15) 



with i,j,l,m = 1,2,3. From the definitions (p],p^ one 
sees that r/ is related to the 12-component of the 2-point 
Green function A^^ of the composite field 7r/„i in the 
single-time representation: 

7?=^ lim All{p)^-l- hm ImAllip). (16) 

ZU p^.p— >0 ZUp'',p^O 

In the second equation we used that the (12)-component 
of any 2-point function is purely imaginary as follows 
from Eqs. (0)-(^. In ( p^ ) each composite field TTim is 
composed of two single particle fields (j). Substituting 
(la) into (nj) one can rewrite (Qq) as 

/3 ,. f d^k f dU , , , 
'?=- „1™„ / 77TZTJ / 77^Z^ Jlm{-k,p+k) 



5pV"o7 (27r)4 7 (27r)4' 

X Im G1122 {-k,p+k, q, -p-q)f"' {q, -p-q) (17) 

where G1122 denotes the (1122)-component of the 4-point 
Green function for the (j) field, Jimip, q) joins two (j) prop- 
agators to a point (see Fig. 1): 



Jlni{p, q) = Piqm - hSlm PQ , 



(18) 



and a symmetry factor 4 accounts for the different pos- 
sibilities of doing so. 

Jeon showed H that in 0* theory at leading order of 
the coupling constant (~ j^-) all planar ladder diagrams 



(see Fig. 1) contribute to the shear viscosity. On the 
other hand, contributions from crossed diagrams can be 
ignored. In the CTP formalism the infinite number of 
ladder diagrams can be resummed by using the Bethe- 
Salpeter (BS) equation. Fig. 2 illustrates the correspon- 
ding BS integral equation for the 4-point Green function, 
with ladders consisting of rungs formed by a one-loop 
diagram connecting two propagators. In general this BS 
integral equation will couple the different components of 
the 4-point Green function in the CPT formalism with 
each other. The key technical issue is therefore whether 
and how these equations can be decoupled. In the fol- 
lowing we will show that in the weak coupling limit this 
is indeed possible, and that the most convenient basis 
for doing so is the r/a representation of the CTP for- 
malism. The previously derived generalized fluctuation- 
dissipation theorem (FDT) ||9|] plays an important role in 
the decoupling procedure. 
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FIG. 1. The planar ladder diagrams contributing to the 
shear viscosity at leading order. 
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FIG. 2. Bethe-Salpeter equation for the 4-point Green 
function. 

From Eq. (||) we deduce 
Gii22{-k,p+k, q, -p-q) = -n{p°) (19) 

X [G1122 {-k,p+k, q, -p-q) + G^aii (-^^ P+^' 9> ^P-q)] 
so that Eq. (|l7| ) can be rewritten as 

V = -^ lim n(pO) / J,™(-fc,p-t-fc) j'"(g, -p-g) 

5 pO,p-.0 Jk^q 

X Im [Gii22+G^2ii] hk,p+k, q, -p-q)] . (20) 

In the physical representation we have (cf. Eq. (M)) 

Im [G1122 + G2211] = —gilm [Grrra + Gaaar (21) 



-rLrrrar~r^aara ^rarr ^araa ^i 
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From the generalized FDT Q we know that only 7 of the 
16 components of Gaia2Q3a4 are independent; choosing 
them as the 4 fully retarded functions plus Garra, Garar 
and Gaarr and using the relations derived in Sec. Ill C 
of Ref. g we find 

Im (G'ii22+G22ii) ~ Im [aGmaa + bGaraa + cGaara 

~rCLL:raaar ~r ^^arra i J^arar ~r Q^aarrj i V j 

where the coefficients a,b, . . . , g are sums and products 
of thermal distribution functions. In the limit p'^^O the 
latter are given by 

a = /3/ Nk{N^-l) + 0((/3p")2) , (23a) 

b = ^-BILLiN',^l)iN^-l) + 0{{(3py) , (23b) 

c = ^ [^.(^,^-1) + A^,(A^^1)] + 0{{f3p'r) , (23c) 

d = ^ [N,{N^~1) - 7V,(7V|-1)] + 0((/3p°)2) , (23d) 

e = ^ m-l) + (A^,'-l)] + 0((/3/)') , (23e) 

/ = ^ [(^^1) + (^1-1)] + 0((/3p")') , (23f) 

g^~Pp°{Ni-l)+0{{f3py), (23g) 

where Nk = l + 2n(fc°). With limpO^o /3p°n(p") == 1 and 
using the symmetry of the integrations over k and q in 
Eq. (|20| ) as well as the symmetry relation (|l^), the shear 
viscosity ( pO| ) is found as 



■^OLlOi^Ct^OiA 






ImG{-k,k,q,-q)ll"\q), 



(24) 



(27r)4 
where h,imik) = -Jim{-k,k) as in ||] and 

(_T — Z ^arra ^aarr • V ^j 

We now set up the BS integral equation for G according 
to Fig. 2. The Feynman rules of the r / a representation 
§ give 

* Gaia2a3ai\ rC, fc, (/, q) ^ V^OJ 

[iA^,aA-k)] [lA^.aAk)] {2n)*5*{k-q) + 

- [iAa,p,{-k)] [iA„.,^,(fc)] (-iA/3i/32/33/34)(-«^7l727374) 



d^s dH 



*A^2&(s)] [*A/33^3(s + ;-fc)] 



(27r)4 (27r)4 

X {i^G/Si^^asaA-iJ^q^-q)] ■ 

Here the factor ^ on the r.h.s. is the symmetry factor 
associated with the bubble connecting the two lines, and 
in the r/a representation the bare 4-point vertex is given 

byi 



4[i-(-iri 



(27) 



Ua is the number of a indices among (ai, a2, a^, a^). 

In general, the full retarded 2-point Green function 
Ara{k) can be expressed as 



Ara{k) 



1 



fc2+ReS(A;) + iImE(fc) 



(28) 



The advanced Green function Aar(k) is given by the 
complex conjugate expression (see (||)). For massless 0^ 
theory ReS(fc) = -XT'^/2A + ©(A^) whereas Imi;(A:) - 
©(A^) |13| , p^ . For weak couphng the propagators thus 
have quasiparticle poles at Rei^^ = ±-\/fc^ + Ar2/24 
with small imaginary parts, and in the limit A — > the 
product Araik)Aar{k) develops a pinch singularity at Ek 
due to the converging pair of poles in the upper and lower 
half of the complex energy plane. Since the products 
Ara{k)Ara{k) and Aarik)Aar{k) havc no pinch singu- 
larities, the corresponding contributions to the shear vis- 
cosity can be neglected for A ^ 1. 

We now evaluate the BS equation for G in this appro- 
ximation. Using (||)-(^ in Eq. ( [2^ ) one obtains for the 
first term on the r.h.s. of (ESh 

Garrai-k, k, q, -q) - ^iV(fc°) A™ (fc) A,, (fc) (29) 



dH 



(27r)4 (27r)' 



Arr{s)Arris+l~k)Gaara{~l, I, 9, ^q) ■ 



Obviously this doesn't decouple: Garra couples to Gaara- 
However, changing the integration variables s — > — s and 
I ^ —I and using the relations (|^) and ( |l^ ) together with 
iV(-fcO) = -iV(fc") one finds that Garrai-k,k,q,~q) is 
an odd function of k. The corresponding contribution 
to the shear viscosity (E4|) thus vanishes by symmetric 
integration. 

For the second term in ( P5| ) we obtain in the same 
approximation from Eq. (Eq) 

Gaarri-k,k,q,-q) = - Ara{k)Aar{k) [l{27:)^SAk-q) 

A2 r -^4. 



d^l 



■Gaarryh h Qj ^Q) 

Aaris)Arais+l-k) 



8 J (27r)4 (27r)4 

Ara{s)Aar{s+l~k) 

+ Arr{s)Arr{s+l — k) 

+Nif){[Ara{s) ~ Aaris)] A„(s+Z-fc) 

+ Arr{s) [Aar{s+l-k) - A™(s+/-fc)]) | 



(30) 



Already here one sees that the BS equation decouples: 
the integral equation involves only the single component 
Gaarr of the 4-point function. To make further progress 
it is convenient to introduce a function M which is ob- 
tained by truncating two of the external legs of the 4- 
point Green function G: 






(31) 



32Q3Q4(, ""J "-J 9l 9) • 

We also introduce the 2-point spectral density 

pik) = i[Araik)-Aarik)] (32) 

and derive from Eq. (H) 

p{k) 



Araik)Aarik) 



2lmY,{k)' 



(33) 



In the pinching pole approximation the first two terms in 
the square brackets in Eq. (BO) can be approximated as 



Arais)Aaris+l-k) + A^r (s) A™ (s+^-fc) 

« [Ara{s) - Aar{s)][Aar{s+l-k) - A™(s+/-fc)] 

= p{s)p{s+k~l). (34) 

With these ingredients, and using Eq. (|^), Eq. (po|) can 
then be simplified as 

Im Mrrrri- k,k,q,-q) = {27rfS^{k-q) + 

dH pjl) pjs) pjs+k-l) ^^ ^g^^ 



4 J (27r)4 (27r)4 ImS(/) 

[l+n(/0)] [l+n(s"+fc"-/")]n(sO) 



1 



i(fc") 



Note that the integral equation remains decoupled when 
expressed through the truncated 4-point function M . 
For the shear viscosity one thus gets 



10 



"^'^ n(fc°)[l+n(fc")]/^(fc) ^^^^ 



(27r)4 
(2^ 



ImE(fc) 
_^/^(q)ImM„rr(-fc,fc, g, -g)- (36) 



This can be brought into the form given by Jeon g by 
defining 



D^ik) 



(2^ 



/^(q)ImMrrrr(-fc,fc,g, -g) ■ (37) 



Inserting Eq. (|3^) into this definition, the BS equation 
derived here is found to coincide with the results obtained 
by Jeon in Eqs. (4.16)-(4.21) of Ref. @. 

Let us summarize: We have shown that a nonpertur- 
bative calculation of the shear viscosity in hot 0"* theo- 
ry can be performed in the real-time (CTP) formalism 
by solving a BS equation. In the weak coupling limit 
the BS integral equation decouples and involves only a 
single component of the real-time thermal 4-point func- 
tion. Our results agree with those of Jeon p|, but our 
derivation is much more compact. This was made possi- 
ble by using the generalized fluctuation-dissipation theo- 
rem which relates the different components of real-time 
thermal n-point functions. 
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